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' Precision measurements, now and at a future linear electron-positron collider (ILC), can provide indirect information 

^Nj , about the possible scale of supersymmetry. Performing a analysis, we illustrate the present-day and possible future 

^jrj' ILC sensitivities within the constrained minimal supersymmetric extension of the Standard Model (CMSSM), varying 

r-{ ^ the parameters so as to obtain the cold dark matter density allowed by WMAP and other cosmological data. The 

current data are in very good agreement with the CMSSM prediction for tan /3 = 10, with a clear preference for 
relatively small values of the universal gaugino mass, ~ 300 GeV. In this case, there would be good prospects for 

, observing supersymmetry directly at both the LHC and the ILC, and some chance already at the Tevatron collider. 

For tan/3 = 50, the quality of the fit is worse, and somewhat larger mi/2 values are favoured. With the prospective 
ILC accuracies the sensitivity to indirect effects of supersymmetry greatly improves. This may provide indirect access 
to supersymmetry even at scales beyond the direct reach of the LHC or the ILC. 
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00 1. INTRODUCTION 

o ■ 

' Measurements at low energies may provide interesting indirect information about the masses of particles that are 
(~| too heavy to be produced directly. A prime example is the use of precision electroweak data from LEP, the SLC, 
the Tevatron and elsewhere to predict (successfully) the mass of the top quark and to provide an indication of the 
O I possible mass of the hypothetical Higgs boson 1] . Predicting the masses of supersymmetric particles is much more 
^ ' difficult than for the top quark or even the Higgs boson, because the renormalizability of the Standard Model and 
• • . the decoupling theorem imply that many low-energy observables are insensitive to heavy sparticles. Nevertheless, 
[ present data on observables such as Mw, sin^ 0cff, (g — 2)^ and BR(6 37) already provide interesting information 
' on the scale of supersymmetry (SUSY), as we discuss in this paper, and have a great potential in view of prospective 
improvements of experimental and theoretical accuracies Q . 

In the future, a linear e'^e~ collider (ILC) will be the best available tool for making many precision measure- 
ments 0. It is important to understand what information ILC measurements may provide about supersymmetry, 
both for the part of the spectrum directly accessible at the LHC or the ILC and for sparticles that would be too heavy 
to be produced directly. Comparing the indirect indications with the direct measurements would be an important 
consistency check on the theoretical framework of supersymmetry. 

Improved and more complete calculations of the supersymmetric contributions to a number of low-energy observ- 
ables such as Mw and sin^ 6'cff have recently become available jj] . These, combined with estimates of the experimental 
accuracies attainable at the ILC and future theoretical uncertainties from unknown higher-order corrections, make 
now an opportune moment to assess the current sensitivities and the projection to the situation at the ILC see 
Ref. for previous studies. In order to reduce the large dimensionality of even the minimal supersymmetric exten- 
sion of the Standard Model, we work here in the framework of the constrained MSSM (CMSSM), in which the soft 
supersymmetry-breaking scalar and gaugino masses are each assumed to be equal at some GUT input scale. In this 
case, only four new parameters are needed: the universal gaugino mass toi/2, the scalar mass mo, the trilinear soft 
supersymmetry-breaking parameter Aq , and the ratio tan /3 of Higgs vacuum expectation values. The pseudoscalar 
Higgs mass Ma and the magnitude of the Higgs mixing parameter can be determined by using the electroweak 
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vacuum conditions, leaving the sign of /i as a residual ambiguity. 

The non-discovery of supersymmetric particles and the Higgs boson at LEP and other present-day colliders imposes 
significant lower bounds on mi/ 2 and toq- An important further constraint is provided by the density of dark matter 
in the Universe, which is tightly constrained by WMAP and other astrophysical and cosmological data 0. These 
have the effect within the CMSSM, assuming that the dark matter consists largely of neutralinos y, of restricting mo 
to very narrow allowed strips for any specific choice of ^0, tan/3 and the sign of /i Thus, the dimensionality of the 
supersymmetric parameter space is further reduced, and one may explore supersymmetric phenomenology along these 
'WMAP strips', as has already been done for the direct detection of supersymmetric particles at the LHC and linear 
colliders of varying energies A full likelihood analysis of the CMSSM planes incorporating uncertainties in the 
cosmological relic density was performed in Ref. |ic| . We extend this analysis to indirect effects of supersymmetry. 

We consider the following observables: the W boson mass, Mvy, the effective weak mixing angle at the Z boson 
resonance, sin^^offi the anomalous magnetic moment of the muon, — (we use the SM prediction based on 
the e+e~ data for the hadronic vacuum polarization contribution, see the discussion in Ref. Q) and the rare h 
decays BR(fe — > 37) and BR(i3s ii^fi^), as well as the mass of the lightest CP-even Higgs boson, Alh, and the 
Higgs branching ratios BR(/i — > bb)/BR{h — > VFVF*). A detailed analysis of the sensitivity of each observable to 
indirect effects of supersymmetry, taking into account the present and prospective future experimental and theoretical 
uncertainties, can be found in Ref. We briefly summarize here the main results on the combined sensitivities at 
present and for the ILC and update some of the results of Ref. 0], which were obtained for rrit ~ 178.0 ± 4.3 GeV, 
using the current experimental central value of rrit — 172.7 ± 2.9 GeV [l]] |. 



2. PRESENT SITUATION 

We first investigate the combined sensitivity of the four low-energy observables for which experimental measure- 
ments exist at present, namely M\y, sin^ ^cft, (ff — 2)^ and BR(6 — s- sj). The branching ratio BR(_Bs for 
which at present only an upper bound exists, has been discussed separately in Ref. 0|. We begin with an analysis 
of the sensitivity to mi/2 moving along the WMAP strips with fixed values of Aq and tan/3. The experimental 
central values, the present experimental errors and theoretical uncertainties are as described in Ref. The exper- 
imental uncertainties, the intrinsic errors from unknown higher-order corrections and the parametric uncertainties 
have been added quadratically, except for BR(6 57), where they have been added linearly. Assuming that the 
four observables are uncorrelated, a fit has been performed with = E„=i (R^n^ - R^^H I is^^f ■ Here K'^'^ 
denotes the experimental central value of the nth observable, so that = 4 for the set of observables included in 
this fit, 1^^"^° is the corresponding CMSSM prediction and ct„ denotes the combined error as specified above. We 
have rejected all points of the CMSSM parameter space with either Mh < 111 GeV [l2j Il3 (taking into account 
theoretical uncertainties from unknown higher orders) or a chargino mass lighter than 103 GeV [l^ . 

The results for tan /3 = 10 are shown in Fig. [T]using mt = 178±4.3 GeV (left) andm* = 172.7±2.9 GeV ^ (right). 
They indicate that, already at the present level of experimental accuracies, the electroweak precision observables 
combined with the WMAP constraint provide a sensitive probe of the CMSSM, yielding interesting information 
about its parameter space. For tan [3 — 10, the CMSSM provides a very good description of the data, resulting in 
a remarkably small minimum value. The fit shows a clear preference for relatively small values of mi/2, with a 
best-fit value of about mi/2 — 300 GeV. This minimum is even more pronounced and located at slightly lower mi/2 
values for the lower top-quark mass value, where especially Mw and sin^ 0off favour lower mi/2 values. The best fit 
is obtained for ^0 0, while positive values of ^0 result in a somewhat lower fit quality. The fit yields an upper 
bound on TO1/2 of about 600 (450) GeV at the 90% C.L. (corresponding to Ax^ < 4.61) for rrit = 178.0 ± 4.3 GeV 
(mt = 172. 7± 2.9 GeV). Some of the principal contributions to the increase in x^ when TOi/2 increases for tan /3 = 10 
are as follows (evaluated for rrit ~ 178.0 ± 4.3 GeV). For Aq = —1TI1/2, mi/2 ~ 900 GeV, we find that {g — 2)^ 
contributes about 5 to Ax^, Mw nearly 1 and sin'^ 6'cff about 0.2, whereas the contribution of BR(6 sj) is 
negligible. On the other hand, for Aq = +2mi/2, which is disfavoured for tan/3 = 10, the minimum in x^ is due 
to a combination of the four observables, but {g — 2)^ again gives the largest contribution for large TOi/2- For 



PSN 0221 



mt = 172.7 ± 2.9 GeV the relative contribution of AIw and sin^ ( 
the Ax^ is more evenly distributed between the observables. 
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to Ax^ increases, so that also for Aq = —mi/2 
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Figure 1: The results of fits based on the current experimental results for the precision observables Mvk, sin^ ^cff , {g — '^)fi 
and BR(6 S'y) are shown as functions of mi/2 in the CMSSM parameter space with CDM constraints for tan/3 = 10 and 
different values of Aq. The left plot shows the results for rrit = 178.0 di 4.3 GeV, and the right plot shows the results for 
mt = 172.7 ±2.9 GeV. 
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Figure 2: The contours in the CMSSM with tan/3 



10 for the mass of the lighter chargino (left), m- + 

^1 



m^o, and the 

^^9 



lighter stau (right), mfj, based on the fits to the parameter space shown in the right plot of Fig. Q (updated from Ref. 01 for 
mt = 172.7 ±2.9 GeV). 



In Fig. 121 the fit results of Fig. ^ (right plot, rrit = 172.7 ± 2.9 GeV) for tan/3 — 10 are expressed in terms of the 
masses of the lighter chargino, ~ ^x°^ ^"^^ ^^"^ mass of the lighter stau, rrif^ (updated from Ref. Q)- The 

best- fit value for is about 200 GeV, while for m^^ the best-fit value is even below 150 GeV. The best-fit values 
for the masses of the neutralinos, charginos and sleptons [il offer good prospects of direct sparticle detection at 
both the ILC ^ and the LHC 15], allowing a detailed determination of their properties [l^]. There are also some 
prospects for detecting the associated production of charginos and neutralinos at the Tevatron collider, via their 
trilepton decay signature, in particular. This is estimated to be sensitive to mi/2 ^ 250 GeV [l^, covering much of 
the region below the best-fit value of mi/2 that we find for tan/3 =10. 

For tan /3 = 50 the overall fit quality is worse than for tan/3 = 10, and the sensitivity to mi/2 from the precision 
observables is lower, as is shown in Fig. O (left). This is related to the fact that, whereas M\y and sin^ 0eff prefer 
small values of rai/2 also for tan/3 = 50, the CMSSM predictions for {g — 2)^ and BR(fe sj) for high tan/3 are in 
better agreement with the data for larger mi/2 values, see the discussion in Ref. 0]. Also in this case the best fit is 
obtained for negative values of Aq, but the preferred values for mi/2 are 200-300 GeV higher than for tan/3 = 10. 
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Figure 3: The results of a fit based on the current experimental results for the precision observables Mw, sin^ Ocb, {g — 2)^ 
and BR(& — » s'y) are shown as functions of mi/2 in the CMSSM parameter space with CDM constraints for different values 
of Ao, nit = 178.0 ± 4.3 GeV and tan/3 — 50 (left). The corresponding contours for the mass of the lighter stau, mf^, are 
shown in the right plot. 



In the case tan/3 = 50 (for rrit = 178 ± 4.3 GeV) the best-fit values for the LSP mass and the lighter stau are still 
below about 250 GeV Q. The fit result for the stau mass is shown in the right plot of Fig. O The minimum for 
the other masses is shifted upwards compared to the case with tan/3 = 10. The best-fit values are obtained in the 
region 400-600 GeV. Correspondingly, these sparticles would be harder to detect. At the ILC with ^/s < 1 TeV, the 
best prospects would be for the production of X1X2 o^' of Tifi. Other particles can only be produced if they turn out 
to be on the light side of the function. 

For the case of tan/3 = 10, all the coloured particles should be accessible at the LHC However, among them, 
only ti has a substantial part of its x^-favoured spectrum below 500 GeV, which would allow its detection at the 
ILC. The same applies for the mass of the A boson. The Tcvatron collider has a sensitivity to TOj^ < 450 GeV [T^ . 
which is about our best-fit value for tan/3 = 10. For the tan f] — 50 case (with nit — 178.0 ± 4.3 GeV) the particles 
are mostly inaccessible at the ILC, though the LHC has good prospects. However, at the 90% C.L. the coloured 
sparticle masses might even exceed 3 TeV, which would render their detection difficult. Concerning the heavy 
Higgs bosons, their masses may well be below ^ 1 TeV. In the case of large tan (3, this might allow their detection 
via the process bb — > bbH/A — > bb t^t" T^. 

While in the fits presented above we kept j4o/mi/2 fixed, we now analyse the combined sensitivity of the precision 
observables M^y, sin^^eff, BR(6 sj) and {g — 2)^ in a scan over the (^1/2,^0) parameter plane. In order to 
perform this scan, we have evaluated the observables for a finite grid in the (Tn-1/2, ^O) 'tt-o) parameter space, fixing 
iTiQ using the WMAP constraint. We restrict ourselves here to the tan/3 — 10 case with mt — 178 ± 4.3 GeV. 
Fig. El (left) shows the WMAP-allowed regions in the (mi/21^0) plane. The current best-fit value obtained via 
fit is indicated.^ The coloured regions around the best-fit value correspond to the 68% and 90% C.L. regions 
(corresponding to Ax^ < 2.30,4.61, respectively). The precision data yield sensitive constraints on the available 
parameter space for mi/2 within the WMAP-allowed region. The precision data are less sensitive to Aq. The 90% 
C.L. region contains all the WMAP-allowed ^0 values in this region of mi/2 values. As expected from the discussion 
above, the best fit is obtained for negative Aq and relatively small values of toi/2. At the 68% C.L., the fit yields an 
upper bound on mi/2 of about 450 GeV. This bound is weakened to about 600 GeV at the 90% C.L. As discussed 
above, the overall fit quality is worse for tan/3 — 50, and the sensitivity to mi/2 is less pronounced (plot not shown 
here). In this case the best fit is obtained for mi/2 ~ 500 GeV and negative ^o- The upper bound on mi/2 increases 



similar analysis within the CMSSM as the one presented here has recently been carried out in Ref. Il9|, where a fit has been 
performed involving all CMSSM parameters. In this analysis, which did not take into account the electroweak precision observables M\y 
and sin^ S^ff , a preference for larger values of tan j3 and somewhat larger mi^2 has been reported. 
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Figure 4: The results of fits for tan/3 — 10 based on the current experimental results (left) and for the anticipated ILC 
precision (right) for Mw, sin^ ^cff, {g ^ 2)^ and BR(6 S7) (including also and BR(/i bb)/'BK{h WW) in the right 
plot) are shown in the (mi/2,^0) planes of the CMSSM with the WMAP constraint. The current best-fit point is indicated in 
both plots. In the right plot two further hypothetical future 'best-fit' values are shown for illustration. The coloured regions 
correspond to the 68% and 90% C.L. regions, respectively. 



3. ILC PRECISION 

We now turn to the analysis of the future sensitivities of the precision observables, based on the prospective 
experimental accuracies at the ILC and the estimates of future theoretical uncertainties discussed in Ref. We 
perform a fit for the combined sensitivity of the observables Mw, sin^^cff, (5 — 2)^, BR(& 37), A-I^ and 
BR{h bb)/BR{h -> WW*) in the (mi/2,Ao) plane of the CMSSM assuming ILC accuracies. The right plot of 
Fig. 2|shows the corresponding fit results for tan/3 = 10. The WMAP-allowed region and the best-fit point according 
to the current situation (see the left plot of Fig. 0J are indicated. Two further hypothetical future 'best-fit' points have 
been chosen for illustration. For all the 'best-fit' points, the assumed central experimental values of the observables 
have been chosen such that they precisely coincide with the 'best-fit' points^. The coloured regions correspond to 
the 68% and 90% C.L. regions around each of the 'best-fit' points according to the ILC accuracies. 

The comparison of Fig. 01 (right) with the result of the current fit (left) shows that the ILC experimental precision 
will lead to a drastic improvement in the sensitivity to rni/2 and Aq when confronting precision data with the 
CMSSM predictions. The comparison of these indirect predictions for mi/2 and Aq with the information from the 
direct detection of supersymmetric particles would provide a stringent test of the CMSSM framework at the loop 
level. A discrepancy could indicate that supersymmetry is realised in a more complicated way than is assumed in 
the CMSSM. For the best-fit values of the current fit the ILC precision would allow one to narrow down the allowed 
CMSSM parameter space to very small regions in the {mi/2, Aq) plane. For the example shown here with best-fit 
values around rni/2 = 300 GeV it is possible to constrain particle masses within about ±10% at the 95% C.L. from the 
comparison of the precision data with the theory predictions. Because of the decoupling property of supersymmetric 
theories, the indirect constraints become weaker for increasing rni/2- The additional hypothetical 'best-fit' points 
shown in Fig. 01 (right) illustrate the indirect sensitivity to the CMSSM parameters in scenarios where the precision 
observables prefer larger values of mi/2- We find that for a 'best-fit' value of mi/2 as large as 1 TeV, which would lie 
close to the LHC limit and beyond the direct-detection reach of the ILC, the precision data would still allow one to 



^We have checked explicitly that assuming future experimental values of the observables with values distributed statistically around 
the present 'best-fit' points with the estimated future errors does not degrade significantly the qualities of the fits. 
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establish an upper bound on mi/2 within the WMAP-aUowed region. Thus, this indirect sensitivity to mi/2 could 
give important hints for supersymmetry searches at higher-energy colliders. For 'best-fit' values of mi/2 in excess of 
1.5 TeV, on the other hand, the indirect effects of heavy sparticles become so small that they are difficult to resolve 
even with ILC accuracies j^l- Whilst the present analysis has been restricted to the CMSSM, similar conclusions are 
expected to apply also in models beyond the CMSSM. Such an analysis is currently in preparation. 
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